Antiferromagnetic order in (Ga,Mn)N nanocrystals: A density functional theory study 
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We investigate the electronic and magnetic properties of (Ga,Mn)N nanocrystals using the den- 
sity functional theory. We study both wurtzite and zinc-blende structures doped with one or two 
substitutional Mn impurities. For a single Mn dopant placed close to surface, the behavior of the 
empty Mn-induced state, hereafter referred to as "Mn hole", is different from bulk (Ga,Mn)N. The 
energy level corresponding to this off-center Mn hole lies within the quantum dot gap near the 
conduction edge. For two Mn dopants, the most stable magnetic configuration is antiferromagnetic, 
and this result was unexpected since (Ga,Mn)N bulk shows ferromagnetism in the ground state. 
The surprising antiferromagnetic alignment of two Mn spins is ascribed also to the holes linked to 
the Mn impurities that approach the surface. Unlike (Ga,Mn)N bulk, these Mn holes in confined 
(Ga,Mn)N nanostructures do not contribute to the ferromagnetic alignment of the two Mn spins. 



I. INTRODUCTION 

Wide band-gap nitride semiconductors are currently 
used in full-color displays, white light sources, and ul- 
traviolet laser diodes for high-density storage systems i 
Such semiconductors combine group-V nitrogen with el- 
ements of group III such as boron, aluminium, gal- 
lium, and indium. Recently, the well-known nitride 
compound GaN has been extensively investigated in the 
form of quantum dots, both with wurtzite^^^ and zinc- 
blende^ crystal structures. The typical phenomena ap- 
pearing in quantum dots are the discretization of the 
electronic spectra and the blue shift of the fundamen- 
tal gaps.^^^ Moreover, the GaN nanocrystals can be 
doped with diluted magnetic impurities such as man- 
ganese. In fact, Ref. 12 shows (Ga,Mn)N quantum dots 
prepared under solvothermal conditions in the wurtzite 
phase. These particles seem to show a ferromagnetic sig- 
nal in the ground state, like bulk (Ga,Mn)N as cal- 
culated for diluted Mn spins. ^^^^ We must note that 
many experiments on bulk (Ga,Mn)N suggest that Mn 
spins are not diluted but forming clusters which give the 
observed ferromagnetism. ^^^^^ A comparative ab-initio 
study between (Ga,Mn)N nanocrystals in wurtzite and 
zinc-blende phases is nevertheless missing. 

In this work we investigate wurtzite and zinc-blende 
GaN quantum dots doped with one or two Mn impu- 
rities within density functional theory. The doping of 
nanocrystals with Mn atoms that replace host cations 
(MuGa) is actually possible, as already confirmed by sev- 
eral experiments on dots larger than ours.^^^^^^^ The ex- 
perimental doping of small nanocrystals of about 1 nm in 
size has not yet been reported. However, it is expected 
from evidence of 2 nm undoped nanoparticles already 
synthesized. Hence, our calculations anticipate future 
experimental work with small doped dots. The theoret- 
ical and computational details are given in Sec. [Ill The 
case of (Ga,Mn)N nanoparticles doped with a single Mn 
impurity is studied in Sec. Illli We show in this section 



that the doping reaction is endothermic and requires high 
temperatures as already confirmed by the experiments.^^' 
The case of (Ga,Mn)N crystallites doped with two Mn 
dopants is investigated in Sec. |IVl where we show that 
the ground-state Mn impurities are antiferromagnetically 
aligned. The antiferromagnetic order of the two Mn spins 
is related to the different chemical environment around 
the empty state induced by the Mn dopant close to sur- 
face. Hereafter, we call "holes" these empty Mn-induced 
states which have mainly Mn d character, but also N.^^^^ 
This hole state lies within the nanocrystal gap near the 
conduction region and therefore does not contribute ef- 
fectively to the ferromagnetic alignment of the two Mn 
impurities. The interesting role of the Mn holes near 
the surface is explained in Sec. [Till and Sec. |IVl Then, 
(Ga,Mn)N quantum dots of about 2 nm in size would be 
likely antiferromagnetic. This possibility is interesting to 
be taken into account for spintronic applications based 
on such III-V Mn-doped nanostructures. Of course, the 
antiferromagnetic behavior of such small dots will also af- 
fect the overall magnetic character of granular solid^~— 
formed by (Ga,Mn)N nanocrystals. 



II. THEORY AND COMPUTATIONS 

We calculate (Ga,Mn)N quantum dots within density 
functional theory, following the Kohn-Sham scheme and 
the projector augmented-wave method, as implemented 
in VASP (Vienna Ab-initio Simulation Package). ^^^^^ 
Apart from the sp valence states inherent in semiconduc- 
tors, we also take into account the Mn states. These 
latter electrons are responsible for the spin-splitted states 
at the gap edges through the sp-d hybridization. For the 
exchange-correlation potential^^ in the Kohn-Sham equa- 
tions we use the generalized-gradient approximation +/7 
(GGA+/7),^^^^ in which U and J are special parameters 
that account for the strong Coulomb and exchange inter- 
actions between the Mn Zd electrons. We take U = A and 
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J = 0.8 as calculated for Mn when doping bulk GaN.^- 

We investigate quasi-spherical nanocrystals of 12 A in 
diameter centered in a Ga position both with wurtzite 
and zinc-blende structures. We passivate the surface dan- 
gling bonds with pseudohydrogens (H*)^^ as an approach 
to quantum dots synthesized in colloidal solutions and 
also grown in semiconductor matrices. These fictitious 
atoms also prevent the appearance of surface states in 
the near-gap spectrum*^ Through passivation, every Ga 
(45^ 4p) dangling bond at the dot surface is attached to a 
pseudohydrogen with a fractional charge of 5e/4, and ev- 
ery N {2s^2p^) dangling bond to a pseudohydrogen with 
a fractional charge of 3e/4. Dot surfaces are perfectly 
saturated and free of defects so that we can avoid any 
perturbation and clearly investigate Mn-Mn magnetic in- 
teractions within nanocrystals. Anyhow, the infiuence of 
surface defects in the magnetic properties of Mn-doped 
dots is an interesting issue for further studies. The doped 
(Ga,Mn)N nanoparticles contain one or two Mn atoms 
which substitute for one or two Ga cations. 

We use the supercell approximation^— to calculate 
wurzite and zinc-blende crystallites which are infinitely 
repeated in space. The size of the supercell is fixed to 
22 A so that surfaces of adjacent dots become separated 
by 10 A and total energies are converged to me Vs. This 
supercell size thus permits an accurate enough descrip- 
tion of Mn-Mn magnetic interactions within nanocrys- 
tals. The atomic positions are fully relaxed until the 
forces on the atoms are small enough (< 0.02 eV/A). The 
input Ga-N bond lengths are taken from bulk GaN in the 
relaxed wurtzite and zinc-blende structures, dca-N = 
1.99 A. For wurtzite GaN, our calculated lattice con- 
stants are a=3.24 A and c=5.29 A; for zinc-blende GaN, 
our lattice constant is a=4.59 A. The cut-off energy in 
the plane- wave basis set is fixed to 500 eV in order to con- 
verge the total energies of bulk GaN and ferromagnetic 
MnN below 1 meV. 




FIG. 1: (Color online) Studied (Ga,Mn)N nanocrystals with 
(a) wurtzite and (b) zinc-blende structures. The Ga atoms are 
drawn with large spheres (pink), the N atoms with medium 
spheres (green), and the H"^ pseudohydrogens with small 
spheres (yellow). The considered positions for the substitu- 
tional Mn impurities in Ga sites are the centers labeled "I", 
the off-center sites labeled "11" , and for zinc-blende dots the 
sites labeled "III". 



III. NANOPARTICLES WITH A SINGLE MN 
IMPURITY 

In this section we study GaN nanocrystals doped with 
a single Mn atom. The dopant concentrations^^ are 
X = 1/13 0.08 for wurtzite quantum dots and x = 
1/19 0.05 for zinc-blende dots. These concentrations 
are in the same order of those used in the experiments.^ - 
The relaxed geometries are depicted in Fig. (U where we 
see that the Mn impurity can occupy the central site 
labeled "I" , the off-center sites labeled "11" , and in zinc- 
blende particles also the sites labeled "III". Positions II 
are near the crystal surface, not exactly at the surface. 
Since results concerning site III are dependent on the 
passivation species, in the following discussions we will 
focus on sites I and II. 



A. Geometric expansion of N atoms around the 
Mn impurity 

As we relax the atomic positions in (Ga,Mn)N quan- 
tum dots, we comment the influence of Mn doping on 
the crystal geometries. In the undoped nanoparticles 
the Ga-N bond lengths calculated for wurtzite and zinc- 
blende structures are 2.02 - 2.04 A for the central Ga 
cations and 1.96 — 2.00 A around the off-center Ga ions 
placed in site II. The smaller bond distances for position 
II as compared with position I are typically due to quan- 
tum dot surfaces. In the doped nanoparticles, the Mn-N 
bonds measure 2.06 — 2.08 A for the central-Mn case and 
1.99 — 2.04 A around the off-center Mn dopants placed in 
site II. Therefore, as compared with the undoped struc- 
tures, the N shell around Mn expands by 2 — 3%, in close 
accordance with the 2% calculated expansion in wurtzite 
and zinc-blende bulk (Ga,Mn)N. We note that this ex- 
pansion was not observed in II- VI (Cd,Mn)Te quantum 
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FIG. 2: (Color online) Substitutional reaction energy for a 
single Mn impurity. Squares (blue) indicate wurtzite struc- 
ture and circles (green) indicate zinc-blende structure. Dotted 
lines are drawn to guide the eye. The most stable position for 
the Mn dopant is site II irrespective of the crystal geometry. 



dots^ and is even contrary to the contraction around 
Mn dopants in III-V (In,Mn)P nanowires.— 



B. Quantum dot stability versus Mn position 

We focus now on the energies involved in the forma- 
tion and doping of (Ga,Mn)N nanocrystals. The cohesive 
energy of the undoped nanoparticle is defined in relation 
to the free atoms, —{Eq — El^)/N, where El^ is the 
energy of the ith free atom, Eq is the total ground-state 
energy, and N is the number of atoms in the dot. Single- 
atom total energies are also obtained within the supercell 
approximation. The computed values are 3.607 eV/atom 
for the wurtzite structure and 3.593 eV/atom for the 
zinc-blende one. In undoped nanocrystals the wurtzite 
phase is thus more stable than the zinc-blende phase. 
However, since the growth of quantum dots (QDs) also 
depends on kinetics and other chemical potentials apart 
from those used in the cohesive energy, the undoped 
GaN crystallites can actually be synthesized both with 
wurtzite and zinc-blende structures. Therefore, we are 
going to study Mn-doped GaN nanoparticles in both kind 
of geometries, wurtzite and zinc blende. 

Anyhow, the energetics of Mn doping for both struc- 
tures is studied by the following reaction. 



GaN QD + Mn 



+2 



(Ga,Mn)N QD + Ga 



+2 



(1) 



in which the Mn+^ dopant substitutes a Ga+^ cation in 
the quantum dot. Reaction (1) points to the main en- 
ergy difference in the process of solvothermal growth dur- 
ing substitutional doping. The cationic form for the Mn 



dopant participates in reactions with the solvent which 
are beyond the scope of this work. For instance, these re- 
actions may include precursors such as GaC13 and MnC12 
in the presence of hexamethyldisilazane (HDMS).^^ To- 
tal energies of ions stem from the calculation of total 
energies of neutral atoms to which the first and sec- 
ond ionization energies (lEs) are summed up, that is, 
£'o(Mn+2)= £;o(Mn)+l^* IE + 2"^ IE, where IE values 
are taken from the literature. The substitutional en- 
ergy A£^i involved in reaction ([1]) is calculated against 
the Mn position and plotted in Fig. [2j Since this en- 
ergy is positive, reaction ([1]) becomes endothermic and 
its activation requires high temperatures. This finding is 
in agreement with the experiments in which (Ga,Mn)N 
nanocrystals are prepared under solvothermal conditions 
at about 350 °C.^^ In addition, the substitutional energy 
A£^i is smaller for site II than for site I. This indicates 
that position II near the surface is more stable than po- 
sition I, in accordance with previous calculations for sub- 
stitutional Mn impurities embedded in II- VI QDs^^ and 
also in III-V nanowires 4^ 



C. Local magnetic moments at the Mn and 
neighbor N sites 

When doping (Ga,Mn)N nanocrystals with a single Mn 
impurity, one 3(i electron of the free Mn atom (3(i^) trans- 
fers to its neighbor N anions and yields a 3(i^ configura- 
tion at the dopant position with a localized Mn hole.^^-i2^ 
Nevertheless, the value of the local magnetic moment at 
the Mn site is different from AjiB due to the sp-d hy- 
bridization. This local moment is 3.99/1^/4. lO/i^ when 
the Mn atom occupies the center of wurtzite/zinc-blende 
quantum dots and 3. 88/1^/3. 88/iB when it is placed off- 
center in position II. Moreover, the Mn dopant induces 
the magnetic polarization of its neighbor N anions. In 
Fig. [3] we show for zinc-blende nanoparticles the spin 
density around the central Mn impurity and its four N 
neighbors. The integration of this density within spheres 
of Wigner-Seitz radii centered in the N atoms results 
in N local magnetic moments of -0.06/iB/-0.09/iB for 
wurtzite/zinc-blende quantum dots. Small modifications 
of these radii would lead to roughly the same local mo- 
ments with the same signs, as already seen for other Mn- 
doped nanocrystals.^^ The exchange coupling between 
Mn and N magnetic moments is thus antiferromagnetic, 
as also calculated for bulk wurtzite and bulk zinc-blende 
(Ga,Mn)N. 

We can compare the Mn local magnetic moments 
obtained for quantum dots and for bulk (Ga,Mn)N 
with similar Mn concentration. For bulk wurtzite 
Gao.917Mno.083N, the computed Mn magnetic moment 
is 3.99/iB; for bulk zinc-blende Gao.937Mno.063N, the Mn 
magnetic moment is 4.03/iB- These bulk values are hence 
similar to the previous Mn magnetic moments centered in 
wurtzite and zinc-blende (Ga,Mn)N nanocrystals. How- 
ever, they show larger differences with respect to those 
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FIG. 3: (Color online) (a) Spin density in zinc-blende quan- 
tum dots doped with a central Mn impurity. The isosurfaces 
correspond to constant densities ±0.05/isA~'^, with light gray 
(red) for the positive value and dark gray (blue) for the neg- 
ative one. (b) The shown spin density is in a plane which 
includes the Mn dopant and two of its four N neighbors. The 
density at the Mn site is chopped for the sake of clarity. See 
how the spherical integration of the spin density around Mn 
and N atoms yields local magnetic moments antiferromagnet- 
ically coupled. 



moments calculated for Mn in site II near the surface. We 
shall see in next Sec. lIIIDl the local densities of states for 
further explanation about such larger differences. 



D. Role of Mn hole in nanostructures 



FIC. 4: (Color online) Local densities of states (LDOS) for 
wurtzite quantum dots doped with a single Mn impurity 
placed (a) in the dot center and (b) in site II off-center. The 
densities of states are projected onto the Mn 3d states in 
black (blue), the neighbor N 2s states with solid black lines, 
and the neighbor N 2p states in gray (orange). The vertical 
dotted lines indicate the Fermi energy. See for Mn off-center 
the hole level placed within the gap at around 2 eV near 
the conduction edge. The zinc-blende densities of states are 
analogous to the wurtzite ones with the only difference of a 
triply degeneracy in the spin-up valence edge state (HOMO) 
for the central-Mn case. Insets widen the conduction edge 
states (ULUMOs). Note the Mn hole level closer to the con- 
duction region in panel (b) as compared with panel (a), and 
the different energy shifts for the up and down ULUMOs. 



We investigate the hole linked to the Mn impurity by 
looking at the local densities of states for Mn in site I 
and site II close to surface. In Fig. |4] we show the local 
projections of the crystal states onto the Mn states 
and N sp states around the Mn dopant. The wurtzite 
case is given in Fig. HJ^a) for a centered Mn impurity 



and in Fig. HJ^b) for an off-center Mn placed in site 11. 
From the densities of states we conclude that (i) the hole 
level associated with the central Mn dopant lies at around 
0.6 eV near the Fermi energy; (ii) nevertheless, the hole 
level associated with the off-center Mn impurity lies at 
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FIG. 5: (Color online) Branch diagram describing the sp- 
d hybridization in a zinc-blende quantum dot doped with a 
centered Mn impurity. On the left we show the splitted Mn 3d 
states in the T^^-symmetry crystal field, in the middle the hy- 
bridized pd bonding and antibonding nanocrystal states, and 
on the right the spin-up 2p states of the four N atoms sur- 
rounding the Mn impurity in the Td-symmetry crystal field. 
The Mn hole is indicated as an open rectangle in the shadowed 
valence region. 



around 2 eV near the conduction region. 

The zinc-blende densities of states are similar to the 
wurtzite densities apart from the central-Mn case which 
shows a triply degeneracy in the spin-up valence edge 
state. Such zinc-blende case is described in detail in 
Fig. m Due to the Td crystal field, the five 3d states 
of the central Mn atom are divided in two groups; one is 
composed of three t2-symmetry states and the other of 
two e-symmetry states. The t2 states of the Mn impu- 
rity hybridize with the p-like t2 states of its neighbor N 
atoms and yield the formation of bonding and antibond- 
ing states. The latter states are degenerate and only 
partially occupied due to the Mn hole. For zinc-blende 
nanoparticles with off-center Mn dopants and also for 
wurtzite dots, the crystal symmetry is T^^-like. 



1. Gap- edge splittings and related sp-d exchange constants 

The sp-d hybridization between the Mn 3d states and 
the sp host states yields an effective Mn-quantum dot 
exchange interaction that splits the crystal states at the 
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FIG. 6: (Color online) The sp-d exchange constants and gap- 
edge splittings for single Mn-doped nanocrystals as a func- 
tion of the impurity position. The open symbols stand for 
Noa exchange constant and the closed symbols stand for No (3. 
Squares (blue) refer to wurtzite structure and circles (green) 
refer to zinc-blende structure. Dotted lines are drawn to guide 
the eye. The bulk exchange constants are indicated on the left 
axis with similar notation. The \Nol3\ values for Mn in the 
dot center (site I) are larger than the bulk ones due to the 
different roles of the Mn holes. 



gap edges. Note that we are interested in the change of 
III-V dot states by Mn impurities. The ab-initio energy 
splittings are explicitly given here and also rewritten in 
terms of sp-d constants, A^o<^ and No/3. Thereby, the 
rescaled splittings are interesting not only for theoreti- 
cians but also for experimentalists, both working on di- 
luted magnetic compounds in the bulk^*"'^^ and in quan- 
tum dots."^^^^ 

The sp-d exchange constants are defined with the fol- 
lowing mean-field theory expressions i^^i^^i^^i^^i^^i^^ 



Noa 



x{Sz] 



and Nop 



(2) 



Here A^o is the number of cations per unit volume 
AE^ = £^^(spin down) — £^^(spin up) is the splitting of 
the up and down Undoped Lowest Unoccupied Molec- 
ular Orbitals (ULUMOs); A£^^ = £;^(spin down) - 
£^^(spin up) is the splitting of the up and down Highest 
Occupied Molecular Orbitals (HOMOs); and (S^) = 4/2 
is the average Mn spin. The splittings AE'^ and AE'^ can 
be extracted from the local densities of states depicted 
in Fig.m The A^o<^ and No/3 exchange constants are cal- 
culated from Eq. (|2]) and presented in Fig. [6] together 
with AE^'^ for different Mn positions. Both wurtzite 
and zinc-blende \No/3\ values are larger for the central 
site I than for site H off-center, as it also happens for 
H-VI quantum dots doped with Mn.^^ 

The dependence of the No/3 exchange constant with 
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the Mn position can be explained in more detail by look- 
ing at the valence states in the gap region. In Fig. [3 
we show these states for two Mn sites, dot center and 
off-center position II. For the central-Mn case, the up- 
down HOMO splitting is large both for wurtzite and 
zinc-blende geometries. However, when Mn is moved off- 
center the spin-up valence levels decrease in energy and 
the spin-down levels slightly increase. These shifts are 
due to a smaller Mn-quantum dot exchange interaction 
caused by a smaller host charge density around position 
II. The up-down HOMO splittings for these Mn sites are 
consequently reduced as compared with those for the dot 
centers, and the N^jS exchange constants become hence 
smaller. We note that for the central-Mn case the charge 
of the up states is mainly distributed over the Mn im- 
purity and its N neighbors, but for the down states it 
suffers a deplection around the Mn dopant that we as- 
cribe to the fact that Mn and N are antiferromagnetically 
coupled. In addition, for Mn in position II the charge in 
the up HOMOs spreads perpendicularly to the x axis and 
globally shows P^-like character. 

The sp-d exchange constants are different for nanopar- 
ticles than for bulk (Ga,Mn)N. In order to assess this 
difference, we calculate N^a and N^jS for two bulk com- 
pounds with Mn concentrations which are similar to 
those in the studied crystallites. For bulk wurtzite 
Gao.917Mno.083N we obtain N^a = —0.23 eV and NqP = 
—5.86 eV; for bulk zinc blende Gao.937Mno.063N we ob- 
tain Noa = 0.28 eV and NqP = -4.47 eV. For compari- 
son, these bulk constants are indicated in Fig. [6] as marks 
on the left axis. The Noa values for Mn in positions I 
and II are comparable to the A^o<^ ones calculated for 
bulk wurtzite and bulk zinc blende. Moreover, the \No/3\ 
values are larger than the bulk ones for Mn in the central 
site I than in site II near the surface. Due to confinement, 
the hole level associated with the central Mn impurity is 
closer to the valence levels within the nanocrystal gap 
than in the bulk. The stronger exchange interaction be- 
tween this Mn hole with d-\ike character and the valence 
edge states increases the splitting of the up and down 
HOMOs and consequently the \No/3\ values. For position 
II as compared with position I, the off-center Mn hole 
lying near the conduction levels indicates a smaller inter- 
action with the valence states, and hence smaller \No/3\ 
values. 



IV. NANOCRYSTALS WITH TWO MN 
IMPURITIES. ANTIFERROMAGNETIC ORDER 
IN THE GROUND STATE 

In this section we investigate (Ga,Mn)N nanoparticles 
doped with two substitutional Mn impurities in the fer- 
romagnetic and antiferromagnetic configurations. The 
ferromagnetic state is calculated for a total magnetic mo- 
ment in the quantum dot (QD) ofSfiB] the antiferromag- 
netic state is calculated for a null magnetic moment in 
the QD. 
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FIG. 7: (Color online) Valence levels in the gap region for 
(a) wurtzite and (b) zinc-blende quantum dots doped with 
a single Mn impurity located in the central site I with open 
rectangles and in an off-center site II with solid lines. The 
charge densities next to the valence states show their global 
S- or P-like character. The geometries with the reference axis 
centered in the dot are oriented so that position II lies on the 
X axis. Crosses indicate site II whenever visible. 
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FIG. 8: (Color online) (a) Substitutional reaction energy for 
two Mn dopants, (b) Interaction energy between two Mn 
impurities. Squares (blue) on the left stand for wurtzite (WZ) 
structure and circles (green) on the right stand for zinc-blende 
(ZB) structure. Closed symbols refer to ferromagnetic Mn 
spins and open symbols refer to antiferromagnetic. Dotted 
lines are drawn to guide the eye. We note from (a) that the 
most stable Mn impurities are aligned antiferromagnetically 
and located in positions I-II both for wurtzite and zinc-blende 
geometries. 



A. Nanocrystal stability versus positions and 
magnetic alignments of the two Mn spins 

By doping with two Mn ions we replace two Ga cations 
as described by the following reaction: 



GaN QD + 2Mn 



+2 



(Ga,Mn)N QD + 2Ga 



+2 



(3) 



The required energy for this double substitution is re- 
ferred to as AE2 and plotted in Fig. [8ja) for different 
Mn-Mn positions and magnetic alignments. As com- 
mented previously, the positive substitutional energies 
indicate that reaction (3) is endothermic and thus ac- 
tivated by increasing the temperature, as it occurs in 
the experiments Fig. [8] also shows that the most sta- 
ble Mn impurities are aligned antiferromagnetically and 
placed in sites I-II close to surface. The calculated an- 
tiferromagnetic ground state was unexpected, since it is 
different from the ferromagnetic alignment of Mn spins 
in bulk (Ga,Mn)NJ^"-i^ We relate it to the different role 
of the hole linked to the Mn dopant that approaches the 
surface. The energy of this hole lies within the nanocrys- 
tal gap near the conduction levels. Therefore, this Mn 
hole does not contribute effectively to the ferromagnetic 
order of the two Mn spins. As a consequence, the Mn-Mn 
coupling becomes antiferromagnetic in the ground state 
as it occurs for two Mn spins in II- VI (Cd,Mn)Te QDs."^^ 
We next study the interaction energy between the two 
Mn dopants which is defined as 



A^int = A^i (Mm ) + A^i (Mn2 ) - AE2; (4) 

Mni stands for the first Mn atom and Mn2 stands for 
the second Mn atom. The calculated interaction ener- 
gies are plotted in Fig. [8l^b) as a function of the positions 
and magnetic couplings of the two Mn spins. The in- 
teraction energy A^int quantifies the relative stability of 
nanoparticles doped with one or two Mn impurities. Pos- 
itive interaction values indicate that the two Mn dopants 
tend to occupy the same nanocrystal and negative values 
indicate that they tend to dope two different quantum 
dots individually. For instance, from Fig. [D^b) it can be 
seen that one crystallite with two antiferromagnetic Mn 
spins in positions I-II is more stable than two nanoparti- 
cles with two single Mn spins placed in sites I and II. We 
stress again that in doping reactions kinetics and other 
chemical potentials different from those used in the co- 
hesive energy could also play an important role and may 
modify previous results concerning the stability of the 
nanostructures. 

Since we are dealing with two Mn impurities, we cal- 
culate their local magnetic moments and total energies 
in the ferromagnetic and antiferromagnetic states. They 
are not explicitely given here but rewritten in terms of an 
effective Mn-Mn exchange interaction, quantified by J^^, 
which is interesting not only for theoreticians^i^i^ but 
also for experimentalists^^^^ working on diluted mag- 
netic semiconductors. The J^^ exchange constant stems 
from the Heisenberg-like Hamiltonian^'' H = — 2 J^^SiS2 
and it is thus defined as 



J' 



dd 



FM 



E^AFM 
^0 



,,FM ,,FM 
/^MniA^Mns 



,,AFM,,AFM ' 



(5) 



where E™ is the total energy of the ferromagnetic state, 
^AFM -g total energy of the antiferromagnetic state, 
and /iMni = is the local magnetic moment at the 
Mni site in Bohr magnetons. The J^^ exchange con- 
stants are calculated and plotted in Fig. [9] as a function 
of the Mn-Mn positions. Negative J^^ values mean anti- 
ferromagnetic alignments between Mn spins and positive 
jdd yai^gg^ ferromagnetic alignments. Fig. [9] shows also 
that in the most stable positions I-II the two Mn atoms 
are antiferromagnetically ordered unlike Mn spins in bulk 
ferromagnetic (Ga,Mn)N. 



B. The sp-d exchange constants and Mn holes 

To look at the modification of the III-V dot states by 
Mn impurities, we now investigate the Noa and NqP ex- 
change constants obtained from the calculated spin split- 
tings at the gap edges and Eq. (j2j) . The computed values 
are shown in Fig. [10] as a function of the positions and 
magnetic couplings of the two Mn spins. The main re- 
sults are the following: 
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FIG. 9: (Color online) The J^^ exchange constants for twc 
Mn impurities as a function of the Mn-Mn positions. Squares 
(blue) denote wurtzite structure and circles (green) denot( 
zinc-blende structure. Dotted lines are drawn to guide th( 
eye. For the most stable positions I-II the negative J^^ val- 
ues indicate that the two Mn spins are unexpectedly orderec 
antiferromagnetically, unlike Mn spins in bulk ferromagnetic 
(Ga,Mn)N. 



(i) |A^o<^| and \NqI3\ values are larger for ferromagnetic 
Mn spins in positions I-II than for farther apart Mn spins 
in positions II-II. This decrease is similar to that observed 
for two Mn spins in II- VI (Cd,Mn)Te quantum dots.^^ 

(ii) N^a values are negative for ferromagnetic Mn 
spins. The hole levels associated with the off-center Mn 
impurities lie at high energies within the nanocrystals 
gaps and push the spin-up ULUMOs above the spin-down 
ULUMOs. 

(iii) Positive A^o<^ values are larger for antiferromag- 
netic Mn spins placed in sites I-II than Ni^a for a doped 
nanoparticle with a central Mn atom. The hole level 
hnked to the spin-down Mn dopant located in site II 
pushes upward in energy the spin-down ULUMO and 
thereby increases the spin splitting of the ULUMOs cor- 
responding to the single-Mn case. 

(iv) Comparing Fig. [6] and Fig. [10] we see that the 
largest \NqI3\ values are for wurtzite and zinc-blende 
quantum dots doped with a central Mn impurity. 

As in the previous discussion of Sec. lIIIDl we now com- 
pare the A^o<^ and NqP exchange constants for nanocrys- 
tals and for bulk (Ga,Mn)N. We calculate the bulk ferro- 
magnetic compounds wurtzite Gao.833Mno.i67N and zinc- 
blende Gao.875Mno.125N. The Mn concentration in bulk 
Gao.833Mno.i67N is similar to x 0.15 in wurtzite quan- 
tum dots; that of bulk Gao.875Mno.125N is similar to 
X 0.11 in zinc-blende dots. For bulk wurtzite we ob- 
tain Noa = -1.06 eV and Nop = -2.51 eV; for bulk zinc 
blende we obtain Noa = 0.26 eV and No/3 = -3.39 eV. 
These exchange values are shown for comparison in 
Fig. [To] as marks on the left axis. For wurtzite nanopar- 
ticles in the ferromagnetic configuration of Mn spins, the 
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FIG. 10: (Color online) The sp-d exchange constants and 
gap-edge splittings for quantum dots doped with two Mn im- 
purities. Squares (blue) on the left refer to wurtzite quan- 
tum dots in the ferromagnetic configuration of Mn spins and 
circles (green) on the right refer to zinc-blende dots. Open 
symbols denote A^o^ exchange constant and closed symbols 
denote No/3. The triangles stand for antiferromagnetic Mn 
spins. Dotted lines are drawn to guide the eye. For compari- 
son, the exchange constants for bulk ferromagnetic (Ga,Mn)N 
are indicated on the two left axis with similar notation. The 
\No/3\ values are larger for ferromagnetic nanocrystals than 
for bulk (Ga,Mn)N. 



\Noa\ values are smaller than in the bulk; on the con- 
trary, for zinc-blende quantum dots in the ferromagnetic 
state, the \Noa\ values are significantly larger than in the 
bulk. For ferromagnetic nanocrystals, the wurtzite and 
zinc-blende \Nol3\ values are both larger than the cor- 
responding \No/3\ bulk constants. These differences can 
also be explained by the different roles played by the Mn 
holes in bulk structures and in quantum dots (as seen in 
previous Sec. IIIIPp . It seems that for nanostructures, 
the situation of the dopant hole within the gap must be 
analyzed in detail in order to understand their basic mag- 
netic properties. 



V. COMPARISON WITH MN-DOPED GAAS 
NANOCRYSTALS 

So far, the antiferromagnetic alignment of Mn in GaN 
quantum dots has been the main result of our discussion: 
Role of Mn hole in GaN nanostructures in Sec. IIII Dl and 
magnetic order of two Mn spins in the last section. In 
fact, the Mn atom close to the crystal surface induces im- 
purity states near the conduction levels. Now, to make 
contact with the antiferromagnetic coupling of such Mn 
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FIG. 11: (Color online) Local densities of states (LDOS) pro- 
jected onto the Mn 3d states for wurtzite and zinc-blende 
quantum dots made of (Ga,Mn)As. The two Mn impuri- 
ties are closely placed in positions I-II and ferromagnetically 
aligned. Note that the two Mn holes approach each other and 
the valence edge unlike in (Ga,Mn)N nanoparticles. 



zinc-blende phase. These exchange energies are roughly 
half the bulk value for Mn atoms sitting in close Ga po- 
sitions, ~200 meV,^^ in agreement with the same de- 
creasing tendency already calculated for other magnetic 
nanostructures such as (Cd,Mn)Te nanocrystals.^^ 

It should be clear from the previous discussion of Mn 
in GaN quantum dots that for ferromagnetic Mn spins 
the so-called Mn hole levels must be close to the GaAs 
valence states. Indeed, the local densities of states pro- 
jected onto the Mn 3d states in Fig. [TT] show that the 
Mn holes in GaAs nanoparticles are close together in the 
nearby of the valence region. These Mn holes behave 
hence as in bulk (Ga,Mn)As^^ and mediate the ferro- 
magnetic alignment of the two Mn spins. 

Looking at the Mn hole levels in its relation to the 
valence (conduction) region, we can see a ferromagnetic 
(antiferromagnetic) behavior of Mn impurities in doped 
nanostructures. The difference in this hole position must 
be found in the smaller bond compression around Mn for 
(Ga,Mn)As quantum dots. 



atoms, it is essential to consider other III-V nanoparti- 
cles. 

Experimentally, (Ga,Mn)As nanocrystals can be cre- 
ated by Mn implantation on GaAs followed by thermal 
treatment, and also by annealing (Ga,Mn)As thin films 
grown by molecular beam epitaxy^^ The self-organized 
nanoclusters are analyzed by x-ray spectroscopy,^^ mi- 
croscopic techniques ^1^^ and also by using SQUID 
magnetometry.^^ Magnetic force microscopy measure- 
ments on (Ga,Mn)As nanoprecipitates show ferromag- 
netic features at room temperature.^^ This experimen- 
tal result and the reported antiferromagnetic behav- 
ior of (Ga,Mn)N quantum dots motivates us to search 
for a change in the magnetic order also in (Ga,Mn)As 
nanocrystals. 

To illustrate this other compound, (Ga,Mn)As, with- 
out taking into consideration all the positions for Mn, 
let us focus on I-II sites, i.e. the case shown to be 
clearly antiferromagnetic (AFM) in (Ga,Mn)N. We cal- 
culate (Ga,Mn)As nanoparticles with two Mn dopants 
which replace two Ga atoms in the close interacting sites 
I-II. We study both wurtzite and zinc-blende geometries 
of about 1 nm in diameter as those given in Fig. [TJ The 
computational details in this case are similar to those 
already explained in Sec. |TT1 

As compared with the undoped structures, in the 
doped ones the As shell around the Mn atom close to 
surface is expanded by 2 — 4%, thus it is more expanded 
than the N shell in (Ga,Mn)N nanoparticles. Moreover, 
the computed quantum dots are ferromagnetic (FM) in 
the ground state, with a FM-AFM exchange interaction 
of 122 meV in the wurtzite phase and 104 meV in the 



VI. CONCLUSIONS 

In summary, we have investigated wurtzite and zinc- 
blende (Ga,Mn)N quantum dots doped with one or two 
substitutional Mn impurities within density functional 
theory. We have obtained that wurtzite and zinc-blende 
structures show similar results as they cannot be distin- 
guished up to second neighbors. Anyhow, there are small 
differences between them commented in the text when 
appropiate. For a single Mn dopant in the dot center, the 
calculated \NoP\ values are larger than in bulk (Ga,Mn)N 
with similar Mn concentration. For two Mn dopants, the 
most stable magnetic state is antiferromagnetic, and this 
was unexpected since bulk (Ga,Mn)N exhibits ferromag- 
netism in the ground state. We ascribe this surprising 
effect in (Ga,Mn)N nanoparticles to the holes linked to 
the Mn impurities placed close to surface. These holes 
do not contribute effectively to the ferromagnetic order 
of the two Mn spins. We show that the ferromagnetic 
behavior of bulk (Ga,Mn)N can be changed by reducing 
the crystal size. 

From the antiferromagnetic result on small (Ga,Mn)N 
nanoparticles, it seems possible that larger dots and other 
nanostructures such as thin films with Mn dopants close 
to surface could also be antiferromagnetic in the ground 
state. Therefore, we hope that our results concerning 
antiferromagnetic (Ga,Mn)N quantum dots will encour- 
age further ab-initio calculations and experiments on Mn 
impurities buried near surfaces of semiconductor nanos- 
tructures. Indeed, recent in-progress work on (Ga,Mn)N 
nanolayers suggests the same antiferromagnetic behavior 
for Mn spins. 
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